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We show first reconstructions of the photon-number distribution obtained with a multi-channel 
fiber-loop detector. Apart from analyzing the statistics of light pulses this device can serve as a 
sophisticated post-selection device for experiments in quantum optics and quantum information. 
We quantify its efficiency by means of the Fisher information and compare it to the efficiency of the 
ideal photodetector. 



A major drawback of common detectors of weak light 
fields is their lack of photon-number resolution. Due 
to the nonlinear character of amplifying process the re- 
sponse of such detectors is not sensitive to the strength 
of the input signal. The only two detection events in 
such a case are "click" and "no click" that correspond 
to the presence or absence of the signal. A device capa- 
ble of photon-number resolution would contribute both 
to fundamental research in quantum optics and to imple- 
mentation of quantum communication and information 
protocols. While such devices were recently indeed con- 
structed [EEHSIHi, they require operation under 
extreme conditions at present and therefore did not be- 
come a common laboratory tool yet. Also, their photon- 
number resolution is still limited to only few photons. 

Another way of circumventing this problem is splitting 
the input pulse using a multiport device followed by an 
array of conventional binary detectors as was proposed 
m 0) 13 • I n the ideal case of many output ports the in- 
put pulse gets perfectly split and each photon is detected 
separately. However reasonable performance of such de- 
vice would require a very large number of beamsplitters 
and detectors which would results in a bulky and costly 
detection device. It has been suggested @ to replace the 
complicated multiport device by a fiber loop and a sin- 
gle photodetector, see Fig. ^ After each round-trip part 
of the incoming pulse gets transmitted to a conventional 
binary detector. This results in a time resolved series of 
detections, each of them corresponding to a different out- 
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FIG. 1: Fiber- loop detector 
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put port of the multiport device. Such a multiple photon 
resolving device has recently been built in our laboratory 
flfjj| . The variable ratio coupler inserted at the entrance 
to the fiber-loop delay line, see Fig. ^ is use d to adjust 
the transmission probabilities of the output channels to a 
certain extent and thus tweak the overall resolving power 
of our instrument. 

The purpose of this communication is twofold. First, 
we will analyze the performance of the fiber-loop detector 
and compare it to the efficiency of the ideal photodetec- 
tion device. Second, we will show how to reconstruct 
the photon statistics of the input pulse from the data 
measured at the fiber-loop detector via the maximum- 
likelihood principle and apply this technique to experi- 
mental data. 

The fiber-loop detector is used to count photons con- 
tained in the input pulse in an indirect way. Therefore it 
is natural to relate its performance to that of the "text- 
book photon counter" — the ideal photodetector of quan- 
tum efficiency fj. The concept of the equivalent efficiency 
is extremely useful in this context. The equivalent effi- 
ciency of the fiber-loop detector is defined to be equal 
to the quantum efficiency of the ideal photodetector that 
gives the same measurement error. This concept was 
used by Braunstein and Nemoto 0] to analyze another 
indirect photon-counting device — the homodyne detec- 
tion scheme. Since their motivation was to use the ho- 
modyne detection for signaling they chose the mutual in- 
formation as their measure of the measurement accuracy. 
Our motivation is different. We are rather interested in 
the analysis of the input light statistics. Let p n be the 
true photon-number distribution (i.e. the diagonal ele- 
ments of the density matrix) of the input pulse. After 
N identically prepared pulses have been detected we can 
make an estimate of the photon number distribution p n . 
Due to the finite measurement time and the presence of 
noise it will usually differ from the true distribution. We 
will take the mean quadratic distance d — ~^2 n {p n — Pn) 2 
between p n and p n as our measure of the estimation er- 
ror. Let us note that for density matrices diagonal in the 
Fock basis this distance equals the Hilbert-Schmidt dis- 
tance between the true and estimated states. Its recipro- 
cal divided by N, I = l/(dN), quantifies the information 



2 




1-T, 




Ti 



1-T, 



FIG. 2: Substituting scheme of the fiber-loop detector with 
three output channels 



about the input pulse gained from the measurement per 
pulse. Naturally, this distance depends not only on the 
measurement device but also on the chosen estimation 
strategy. Let us calculate the error of the optimal esti- 
mation: It is well known that no matter how clever esti- 
mation procedure is adopted, its error cannot be smaller 
than the Cramer- Rao bound [H,[l3|]. This bound applied 
to I implies that 
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Now we will analyze the fiber-loop detector which can be 
thought of as a multiport device. Such a multiport, for 
simplicity limited to only three output channels, is shown 
in Fig. H 

The simplest situation occurs when the coincidences 
between detections at different output channels are not 
registered. Then an s-channel device reduces to s in- 
dependent binary detectors of generally different quan- 
tum efficiencies. As was shown in such a detec- 
tion scheme already yields enough information for the 
reconstruction. In this case the independent events are 
zero-detection events in the individual channels having 
probabilities 
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where F is the Fisher information matrix, 



A significant gain in accuracy of the reconstruction can 
be expected when all possible coincidences are taken into 
account. A detection event of such a complex observation 
is best recorded in the binary notation where "1" or "0" 
appearing at a certain position corresponds to a "click" or 
"no click," respectively, in the given output channel. The 
most simple nontrivial case of such a measurement is the 
fiber-loop detector with only two output channels (the 
pulse goes round the loop only once). This is equivalent 
to letting the transmissivity Ti of the second beamsplitter 
in Fig. [21 go to unity. In this case we have four detection 
events. Their probabilities read 
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and P(n\p) is the conditional probability of detecting the 
outcome n provided the true statistics of the input pulse 
is p. The averaging is done over the all possible out- 
comes n of the measurement. The probability P(n|p) 
completely describes the measurement apparatus. Pro- 
vided the individual detection events are independent, 
P(n|p) is a multinomial distribution, 
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Here p n is the probability that the input pulse gives rise 
to a detection in a particular channel n. The Fisher in- 
formation matrix then simplifies to 
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In the case of the ideal detector of efficiency rj the inde- 
pendent outcomes are just the different counted numbers 
of photons and the corresponding probabilities are given 
by the well-known Bernouli distribution, 
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Poo = Y,P^ 1 -inT-v2(i-T)] m , 

m 

Pio = Ep m [l-r? 2 (l-T 1 )] m -Poo, 

m 

Pol = y^Pmjl -m T ) m - POO, 
m 

Pn = 1 - Poo -Pio ~Poi, 



and T = T\. It is clear that only weak input signals can 
fully be analyzed by the simple measurement 10. In the 
following we will assume that the input density matrix 
can be truncated and thus fully specified by giving its first 
few elements. The most simple case is a pulse containing 
at most two photons, 

p=(l- Pl - p 2 )\0)(0\ + pi|l)(l| + p 2 |2)(2|. (8) 

A successful estimation consist in identifying the single 
photon and two photon contributions. Using Eqs. fijl. 

and Q the information about the state (jHJ 
yielded by the fiber-loop detector can easily be calcu- 
lated. As the resulting expression is rather complicated 
and not suitable for discussion we will consider only the 
leading term of its expansions in rjs. Even best binary 
detectors have rj significantly smaller than unity guar- 
anteeing the rapid convergence of the series. Setting 



FIG. 3: Performance of the fiber-loop detector with two chan- 
nels for various input states and splitting ratios. Curves from 
below: p x = 0.9, p 2 = 0.1; pi = 0.98, p 2 = 0.2; p x = 0.998, 
p 2 = 0.002. 



FIG. 4: Equivalent ideal detection efficiency of the fiber-loop 
detector with three output channels in dependence on the 
quantum efficiency of the binary detector used. The efficiency 
of standard commercial single-photon detectors is about 60% 
(dashed line). 



Vi = ? 72 — V the information reads, 
2r?7 2 

/loop ~ -= — , P2 # 0, (9) 

/loop » — + 2(1-t)t7 2 , p 2 = 0, (10) 
Pi 

and t = T(l — T). As can be seen the information is de- 
termined by the last non- vanishing p n . Naturally, states 
with large multi-photon contribution are more difficult 
to estimate. If p% significantly differs from zero it is opti- 
mal to split the input pulse between the output channels 
equally, T = 1/2, see also Fig. (J3J. However for states 
with negligible two-photon content p 2 ~ the informa- 
tion exhibits a local minimum at this point, see Fig. 
This minimum is caused by the term in Eq. (|10f> which is 
quadratic in rj and since rj is always significantly smaller 
than unity the minimum s a shallow one. In an experi- 
ment where the true state is not known, T = 1/2 would 
be the best choice. The corresponding information gain 
of the ideal detector is 

fj 2 

/ideal ~ ~ , P2 ^ 0, (11) 

5^2 

/ideal PS — , P 2=0. (12) 
Pi 

As can be seen, for states with negligible two-photon con- 
tent the ideal and fiber-loop detectors show the same per- 
formance. This is not surprising because for single pho- 
ton states there is no difference between the ideal and 
binary detector. However for states with significant two- 
photon content the information the fiber loop detector 
becomes inferior to the ideal one. The equivalent effi- 
ciency of the fiber-loop detector yields in the worst case 
can be calculated by comparing Eqs. © and (|ll|l . 

r hq « 77/V2. (13) 

The dependence of the equivalent efficiency calculated by 
the exact inversion of the Fisher matrix on the quantum 



efficiency of the binary detector used in the fiber-loop 
detector is shown in Fig. As one can see Eq. 113|) 
makes an excellent approximation. Let us close this part 
noting that for weak input fields the information gain 
of the fiber-loop and ideal detectors differ only by the 
factor of v2< Further improvement can be achieved by 
increasing the number of output channels. Provided the 
noise of the detectors can be neglected, splitting output 
channels into two or more new channels leads to a refine- 
ment of the probability operator measure describing the 
detector. Any such refinement can only increase the in- 
formation gain. Such many-channel fiber-loop detectors 
will be discussed elsewhere. 

In real experiment experimenter's control of the pa- 
rameters of the fiber-loop detector is severely limited for 
many reasons, so it is not possible to choose the optimal 
configuration of the device. 

Prior to its utilization the device has to be calibrated. 
We use coherent pulses whose statistics is known to be 
Poissonian for this purpose. From the probabilities of 
zero counts at the individual output channels the trans- 
missions of beamsplitters in the substituting scheme are 
readily calculated. In the next step we calculate proba- 
bilities p n of all possible outcomes of the experiment. We 
start from the probabilities of "zero-arbitrary" detection 
events. Such events consist of detecting nothing in given 
output channels while we do not care about the rest; they 
read 

dim s 
fc=0 1=0 

where dim is the chosen cutoff of the Hilbert space and 
€ji is zero if there was an arbitrary event (denoted V) 
detected at the l-th channel and unity if there was no 
detection, and r\\ = (ni)/N is the detection efficiency 
of the l-th channel. Now the probability of having a 
coincidence at the first and third channels of a three- 
channel detector can be written in terms of the single- 
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FIG. 5: (a) Reconstruction of the photon-number distribu- 
tion. Black dots are the best Poissonian fit with (n) ~ 0.6 
photons per pulse, (b) The calculated equivalent efficiency of 
the fiber-loop detector used in this experiment. 
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count probabilities as follows: pioi = Pvov — Pioo — Pool — 
Pooo- The latter can, in turn, be expressed as pioo = 
Pvoo - Pooo and p 001 = p oov - p 000 , respectively. The 
probabilities of all 2 s detection events possible with an 
s-channel loop detector can be easily calculated using 
this simple recursive procedure. Having p n determined 
one can proceed with the reconstruction. Notice that p n 
are linear combinations of density matrix elements p n , 

Pj =^2cjiPi, j = 1, . . . ,2 s , (15) 

i 

Cji = dpj/dpi being defined by the parameters of the 
fiber-loop detector. The problem l|15|) is a typical lin- 
ear and positive problem. The latter property follows 
from the non- negativity of p n . Of course, the experi- 
mentally observed frequencies will generally differ from 
the theoretical probabilities, fj ^ pj. Because of this 
the problem 115|) has usually no exact solution and has 
to be solved in the statistical sense, for example via 
minimization of some statistical measure of distance be- 
tween fj and pj. We use the approach based on mini- 
mizing the Kullback-Leibler divergence (also called rel- 
ative entropy) d = T,jfj lo EPj ~ T,jfj lo Efj between 
data and theory. The state p n minimizing d is at the 
same time the maximum-likelihood estimator of the in- 
put state and therefore attains the Cramer-Rao bound 
asymptotically for large N . The minimization can be car- 
ried on e. g. v ia the iterative Expectation-Maximization 
algorithm |15lllq: 

p^)=p(»)^^M_. (16) 

Left panel in Fig. [5] shows the reconstruction of the 
photon statistics of a laser pulse containing 0.6 photons 
on the average. The agreement between the measured 
and theoretical statistics is very good. In this case only 
first three output channels were sufficient for the recon- 
struction. The equivalent efficiency of the three-channel 
detector calculated for Poissonian light is shown on the 
right in Fig. |SJ We would like to emphasize that this is 



FIG. 6: Reconstructed photon- number distributions of dif- 
ferent sources of light pulses (bars); best theoretical fits 
(black dots); in panel (c) best Poissonian fit is shown for 
comparison (triangles). Fluctuations: (a) stable amplitude, 
P(I) oc 5(1 — To); (b) two-component amplitude, P(I) oc 
5(1— Ii)-{-5(I— 12); (c) uniform over the range of [I3 — 5, 13 + 5]; 
(d) exponential, P(I) = exp(— I/I^j/I^. Resulting statistics: 
(a) Poissonian; (b) composite Poissonian; (c) regularized gen- 
eralized incomplete Euler gamma function; (d) Bose- Einstein. 
Fitted parameters: I = 4.6, h = 0.94, h = 4.6, I 3 = 2.0, 
I4 — 2.1, and 5 — 1.2 photons per pulse. 

the efficiency of the real experimental apparatus where a 
significant fraction of light is either lost inside the appa- 
ratus or not detected due to small quantum efficiency of 
detectors Nonetheless, in a wide range of intensities 
our detector is the equal of the ideal detector of quantum 
efficiency of about 25%. 

More intense pulses require the use of a larger num- 
ber of output channels. We have tested our fiber-loop 
detector on different kinds of light pulses. The statistics 
has been artificially changed by changing the intensity 
of subsequent Poissonian pulses using a random num- 
ber generator. The statistics P(I) of intensity fluctua- 
tions determines the true photon statistics of light pulses 
through the composition rule 

/jn 
—e-'P^dl. (17) 

Experimental results obtained with different light sources 
are summarized in Fig. Notice that the composite na- 
ture of p nicely shows in the upper right panel. Also the 
super-Poissonian character of the input pulse is clearly 
seen in the lower left panel. Several tens of thousands of 
pulses were used for each reconstruction. 

In conclusion we have demonstrated the photon- 
counting capability of the fiber-loop detector and quan- 
tified its efficiency by means of the Fisher information. 
The possibility to reconstruct the statistics of the input 
light pulses of moderate intensities have also been shown. 
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